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Abstract Recent observations as well as theoretical studies of YSO jets suggest the 
presence of two steady components: a disk wind type outflow needed to explain the 
observed high mass loss rates and a stellar wind type outflow probably accounting 
for the observed stellar spin down. In this framework, we construct numerical two- 
component jet models by properly mixing an analytical disk wind solution with a 
complementary analytically derived stellar outflow. Their combination is controlled 
by both spatial and temporal parameters, in order to address different physical con- 
ditions and time variable features. We study the temporal evolution and the interac- 
tion of the two jet components on both small and large scales. The simulations reach 
steady state configurations close to the initial solutions. Although time variability is 
not found to considerably affect the dynamics, flow fluctuations generate condensa- 
tions, whose large scale structures have a strong resemblance to observed YSO jet 
knots. 



1 Introduction 

In the last few years, a promising two-component jet scenario seems to emerge in 
order to explain Young Stellar Object (YSO) outflows. Observational data of Clas- 
sical T Tauri Stars (CTTS) [1], [3] indicate the presence of two genres of winds: one 
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being ejected radially with respect to the central object and the other being launched 
at a constant angle with respect to the equatorial plane (e.g. Tzeferacos et al., this 
volume). In turn, CTTS outflows may be associated with either a stellar origin, or 
a disk one or with both wind components having roughly equivalent contributions. 
In addition, such a scenario is supported by theoretical arguments as well (e.g. [2]). 
An extended disk wind is required for the explanation of the observed YSO mass 
loss rates, whereas a pressure driven stellar outflow is expected to propagate in the 
central region, being a strong candidate to address the protostellar spin down [6]. 

The goal of the present work is to study the two-component jet scenario, taking 
advantage of both analytical and numerical approaches. In particular, we construct 
numerical models by setting as initial conditions a mixture of two analytical YSO 
outflow solutions (each one describing a disk or a stellar jet), ensuring the dom- 
inance of the stellar component in the inner regions and of the disk wind in the 
outer. The combination is achieved with the introduction of few normalization and 
mixing parameters, along with enforced time variability of the stellar component. 
We investigate the evolutionary properties, steady states and the features of the final 
configurations of the dual component jets. Although the detailed launching mecha- 
nisms of each component are not taken into account, such models seem capable to 
capture the dynamics and describe a variety of interesting scenarios. 

The employed analytically derived MHD outflows, defined as ADO (Analytical 
Disk Outflow; denoted with subscript D) and ASO (Analytical Stellar Outflow; de- 
noted with subscript S), have been derived in the context of self-similarity [9] and 
each one effectively describes a disk wind [10] or a stellar jet [8], respectively. In 
Matsakos et al. [4], we have addressed the topological stability, as well as several 
physical and numerical properties, separately for each solution. This article summa- 
rizes the numerical setup and reports the results of few significant cases of the dual 
component jet. A thorough study can be found in Matsakos et al. [5]. 



2 Numerical two-component jet models 

The two-component jet model parameters can be classified in two categories. The 
first one contains those associated to the relative normalization of the analytical 
solutions, i.e. the ratios of the characteristic scales of each model (denoted with 
subscript *, calculated on a specific fieldline at the Alfvenic surface): 

4 = ^, ty = ^, 4 = |^, (1) 
r* Vd* i>D* 

where is the Alfvenic spherical radius of the ASO model, r* is the cylindrical 
radius of the Alfvenic surface of the ADO model (of a specific fieldline) and the 
subscripts L, V and B stand for length, velocity and magnetic field, respectively. 
We assume that the protostar has a solar mass and a radius of ~ 0.01 AU. Since the 
disk wind launching region lies in the range 0.2 — 3AU, we derive £l = 0.1 and 
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ly = 5.96. On the other hand, we define £g = 2, which is the parameter controlling 
the relative dominance of each model. 

The second class of parameters concerns the mixing. In particular, we choose the 
combination to depend on the magnetic flux function A, which labels the fieldlines 
of each analytical model (Ap or A$). Therefore, we define a common trial magnetic 
flux A tr =Ad+As and then all physical variables U are initialized with the help of 
the following mixing function: 



Ujcomp = < 1 - exp 



Atr 
qA„ 



■ t/o + exp 



Atr_ 

qA m 



Us, 



(2) 



where A m is a constant corresponding to the matching surface rooted at 0.16 AU, q 
is a parameter that effectively moves this surface closer to the protostar and d sets 
the steepness of the transition from the inner ASO to the outer ADO solution. We 
choose q = 0.2 and d = 2, whereas a complete parameter study (including £g) can 
be found in [5]. Essentially, Eq. 2 provides an exponential damping of each solution 
around a particular fieldline of the combined magnetic field. 

Moreover, since accretion and protostellar variability are expected to introduce 
fluctuations we multiply the inflow velocity with the following function: 



f s (r,t) = 1 + ^sin 



1 var J 



exp 



r 

2r„ 



(3) 



where T var is the period of the pulsation and r m is roughly the cylindrical radius at 
which the matching surface intersects the lower boundary of the computational box. 
Outflow variability produces the formation of knot-like structures: the introduction 
of radiation cooling (Tesileanu et al. this volume) will allow direct comparison with 
observational data. 

The simulations are performed with PLUTO 1 (Mignone et al. [7]). A uniform 
resolution of 256 zones for every 100 AU is used whereas the simulations have been 
carried out up to a final time of 80 y. On the lower boundary we keep fixed all vari- 
ables to their initial values, on the axis we apply axisymmetric boundary conditions 
and at the upper and right borders of the domain we prescribe outflow conditions. 



3 Results 



In the left panel of Fig. 1, the initial setup (left) and the final configuration (middle) 
of the two-component jet are displayed. The model shows remarkable stability and 
reaches a well defined steady state in only a few years. In particular, the disk wind 
solution remains almost unmodified whereas the stellar component gets compressed 



A versatile shock-capturing numerical code suitable for the solution of high-Mach number flows. 
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Fig. 1 Left panel: Logarithmic density contours (the code unit is 10~ I2 gcm~ 3 ) and magnetic 
fieldlines for the initial two-component jet model (left), final steady state (middle) and when a 
monthly flow variability is applied (right). Right panel: The quantity I0 3 p 2 \/T (roughly related 
with emissivity) is plotted for the yearly variable model. Although max(10 3 p 2 \/r) = 53.9, the 
color bar uses a lower maximum value to enhance the displayed features. 

around the axis. Moreover, a shock manifests during time evolution, located roughly 
along the diagonal line which crosses (10, 40) and (30, 100) (steady state plot). 
This shock is found to causally disconnect the acceleration regions from the jet 
propagation physics and the subsequent interaction with the outer medium. Note 
that there is no such "horizon" present in the initial setup. Furthermore, on the right 
plot of the left panel of Fig. 1, the same model is displayed when a monthly time 
variable velocity is applied on the lower boundary (Eq. 2). Evidently, despite the 
strong gradients seen in the density and the wiggling of the magnetic fieldlines, the 
general structure is retained, proving the stability of the two-component jet model. 

The fact that the system remains very close to the initial configuration, demon- 
strates that the analytical solutions provide solid foundations for realistic two- 
component jet scenarios. Consequently, specific YSO systems can be addressed 
more accurately by constructing analytical outflow solutions with the desirable char- 
acteristics, before merging them into a two-component regime. 

On the right panel of Fig. 1 a quantity related to emissivity is plotted in larger 
scales when a yearly variability is applied. Near the base, the numerical solution 
remains close to the initial ADO and ASO models. However, farther away along 
the flow the fluctuations create knot-like structures, which may be related with jet 
variability. In fact, note that the model is associated with a condensations spacing 
~ 100 AU, similar to the knot spacing of HH30. 

Finally, although not presented in this article, an other important parameter is the 
one controlling the relative contribution of each component, Ib, with which we can 
effectively and smoothly switch the model from a totally magneto-centrifugal wind 
to a pressure driven jet [5]. 
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To sum up (taking also into account the results of [4] and [5]), most of the techni- 
cal part concerning two-component jets, e.g. 2.5D stability, steady states, parameter 
study, time variability etc., is now at hand [5], providing us with all the necessary 
ingredients to address YSO jets. Namely, with a) the proper analytical solutions, i.e. 
desirable lever arm, mass loss rate etc., b) the correct choice of the mixing param- 
eters and c) an enforced time variability that effectively produces knot structures, 
we are now ready to qualitatively study different and realistic scenarios, address 
observed jet properties and ultimately understand the various outflow phases of spe- 
cific T Tauri stars. 
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